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RESULTS  FROM  A  TWO  LAYER  AND  A  REDUCED 
GRAVITY  MODEL  OF  THE  ALBORAN  SEA 

INTRODUCTION 

The  circulation  of  the  Alboran  Sea  Is  dependent  on  both 
Atlantic  and  Mediterranean  waters.  Atlantic  water  flows  through  the 
narrow  (20  km  wide)  and  shallow  (300  m  deep)  Strait  of  Gibraltar  Into 
the  Alboran  Sea  forming  a  150-200  m  deep  surface  layer  (Ovchinnikov, 
1966;  Lanolx,  1974;  Katz,  1972).  Medl terranean  water  enters  the  Alboran 
at  Its  open  eastern  boundary  and  flows  slowly  westward  In  the  form  of  an 
Intermediate  and  lower  layer.  It  has  been  suggested  that  even  the 
deepest  water  can  exit  through  the  Strait  (Stommel  et  al.,  1973).  The 
large  volume  transport  of  Inflowing  Atlantic  water,  1  to  2  X  106  m3/sec, 
(Lacombe,  1971;  Bethoux,  1979;  Lacombe,  1982)  retains  Its  Identity  as  a 
narrow  (30  km  wide)  jet  with  Initial  speeds  of  100  cm/sec  near  the 
Strait  (Lacombe,  1961;  Peluchon  and  Donguy,  1962;  Grousson  and  Faroux, 
1963;  Lanolx,  1974;  Cheney,  1977;  Petit  et  al.,  1978;  and  Wannamaker, 
1979).  The  jet  enters  the  basin  and  flows  northeast  to  approxlmatley 
4*W,  curves  southward  and  then  splits  (Figure  1).  Part  of  the  jet  flows 
to  the  west  and  Is  incorporated  In  an  antlcyclonlc  gyre,  while  the 
remainder  flows  southeast  to  Cape  Tres  Forcas  and  then  along  the  African 
coast  forming  the  southern  periphery  of  a  cylconlc  circulation. 

Satellite  Infrared  Imagery  Indicates  variations  In  the  shape, 
location  and  Intensity  of  the  persistent  antlcylconlc  gyre  which 
dominates  the  western  Alboran  Sea.  Hydrographic  data  and  satellite 
Infrared  Images  support  the  year-round  persistence  of  the  gyre,  although 
Its  size  and  location  varies  (Stevenson,  1977;  Cheney,  1978;  Wannamaker, 
1979;  Burkov  et  al.,  1979;  Gallagher  et  al.,  1981).  In  the  eastern 
portion  of  the  Alboran  Sea  a  general  pattern  of  alternating  cyclonic  and 
antlcyclonlc  circulations  has  been  observed  (Cheney,  1978;  Lanolx, 
1974).  Satellite  Imagery  shows  that,  compared  to  the  western  Alboran, 
this  circulation  pattern  Is  far  more  variable  and  of  smaller  scale. 


Figure  1.  Dynamic  topography  of  the  surface  relative  to  200  dbar  for  July-August  1962 
Lanov  (1974). 


A  study  of  circulation  In  the  Alboran  Sea  begins  by  using  the 
simplest  model  capable  of  simulating  major  features  of  the 
circulation.  In  a  system  with  two  so  clearly  defined  layers,  a  reduced 

gravity  model  may  be  used  to  determine  the  circulation  of  the  upper 

layer,  and  a  two  layer  model  to  determine  the  entire  system.  The 

reduced  gravity  model  allows  barocllnlc  instability  and  the  Inclusion  of 
bottom  topography.  Previous  results  using  a  semi-implicit  reduced 

gravity  model  with  a  rectangular  geometry  are  discussed  In  Preller  and 
Hurlburt,  (1982).  We  have  now  progressed  to  a  reduced  gravity  model 
with  real  geometry  and  a  two  layer  model  with  geometry  and  bottom 
topography . 
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THE  MODEL 


The  non-linear  reduced  gravity  and  two  layer  models,  developed 
for  the  Gulf  of  Mexico  by  Hurlburt  and  Thompson,  (1980),  have  been 
adapted  and  modified  for  the  Alboran  Sea.  The  reduced  gravity  model  Is 
an  explicit  variant  of  the  Hurlburt  and  Thompson  model  and  Is  described 
In  detail  by  Preller  and  Heburn,  (1983).  The  two  layer  model  Is  solved 
numerically  using  the  economic  semi -Implicit  method.  The  reduced 
gravity  model  consists  of  an  active  upper  layer  and  a  lower  layer  which 
Is  Infinitely  deep  and  at  rest.  The  two  layer  model  has  an  upper  layer, 
a  finite  lower  layer  and  may  Include  bottom  topography  (Figure  2).  Both 
models  are  stably  stratified  and  have  a  fixed  density  contrast  between 
two  Immiscible  layers.  These  models  assume  a  hydrostatic,  Bousslnesq 
fluid  In  a  rotating  right-handed  coordinate  system  on  a  -plane.  The 
vertically  Integrated  model  equations  are: 
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i  ■  1,  2  Indicates  upper  and  lower  layers 

and  x  and  y  are  tangent-plane  Cartesian  coordinates  with  x  directed 
eastward  and  y  northward,  and  v^  are  the  eastward  and  northward 
velocity  components  in  each  layer,  h^  Is  the  thickness  of  each  layer,  t 
is  time,  g  is  the  acceleration  due  to  gravity,  Is  the  density  of 
seawater  In  layer  1,  fQ  and  yQ  are  the  values  of  the  Coriolis  parameter 
and  the  y-coordlnate  at  the  southern  boundary,  t.  Is  the  wind  stress, 
and  t2  Is  the  Interfacial  stress. 

Figure  3  depicts  the  Irregular  geometry  and  topography  used  In 
the  two  layer  experiments  and  the  explicit  reduced  gravity  cases.  This 
Irregular  boundary  was  created  from  the  Synthetic  Bathymetric  Profiling 
System  (SYNBAPS)  (Vanwyckhbuse,  1973  and  Vanwyckhouse,  1979)  data  set  of 
topography.  These  values  were  Interpolated  from  the  10  minute 

resolution  of  the  model.  The  coastline  was  then  drawn  using  a  cutoff 
depth  of  20  m  to  eliminate  errors  existing  on  the  boundary  of  the  data 
set.  The  four  point  Island  In  the  center  of  the  basin  Is  an  Idealized 
representation  of  the  Alboran  Island.  A  four  point  Island,  20  km  X  10 
km.  Is  the  smallest  allowed  by  the  numerical  schemes  used  In  the 
model.  Except  at  the  ports,  boundaries  are  rigid  and  a  no  slip  boundary 
condition  Is  prescribed  on  the  tangential  flow.  Velocity  Is  prescribed 
at  the  inflow  port.  Normal  flow  at  the  outflow  port  Is  self-determined 
with  the  Integral  constraint  that  the  total  transport  out  of  each  layer 
match  the  Inflow.  In  the  two  layer  cases,  the  amplitude  of  the 
tangential  component  was  set  to  zero  one-half  the  grid  distance  outside 
the  physical  domain  at  both  ports.  A  variation  of  the  Orlanskl  boundary 
condition  (see  Preller  and  Heburn,  1983)  was  used  on  outflow  In  the 
explicit  reduced  gravity  model. 
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Model  parameters  for  each  of  the  model's  standard  cases  are 
given  in  Table  1.  The  only  dissipation  in  the  model,  other  than  the 
viscous  boundary  layer,  Is  due  to  horizontal  friction.  The  value  for 
the  coefficient  of  eddy  viscosity  was  chosen  to  keep  the  grid  interval 
Reynol ds  number 

Re  ■  IjL  30, 

thus  keeping  the  effects  of  frictional  damping  from  dominating  the 
solutions. 

The  value  chosen  for  reduced  gravity  due  to  stratification  was 
appropriate  for  observed  density  differences  in  the  Alboran  Sea.  In  the 
two  layer  model,  the  value  of  g'  In  the  table  Is  multiplied  by 
(Hj  +  H2)/H2  to  yield  the  same  Internal  values  for  gravity  wave  speed, 
radius  of  deformation  and  Rossby  wave  speed  as  In  the  reduced  gravity 
model.  The  undisturbed  upper  layer  thickness  Is  based  on  many 
observations  of  the  Alboran  Sea  which  show  a  sharp  density-salinity 
gradient  at  150-200  m  depth. 

Horizontal  grid  spacing  was  chosen  to  allow  the  resolution  of 
the  mesoscale  features  being  studied  In  the  Alboran  Sea.  A  minimum  of 
8-10  grid  points  Is  necessary  to  properly  resolve  a  feature.  The 
Alboran  Gyre  Is  defined  with  an  average  diameter  of  100  km.  Thus  the 
model  has  a  minimum  of  10  grid  points  with  which  to  define  this  feature. 

For  cases  using  the  Irregular  geometry,  the  western  Inflow  port 
Is  20  km  wide,  an  appropriate  width  for  the  Strait  of  Gibraltar,  while 
the  eastern  outflow  port  Is  245  km  wide.  In  the  lower  layer  of  the  two 
layer  cases,  the  eastern  boundary  serves  as  the  Inflow  port  while  the 
Strait  of  Gibraltar  becomes  the  outflow  boundary.  Inflow  is  exactly 
compensated  by  outflow  separately  in  each  layer.  Specification  of  the 
model  forcing  Is  accomplished  using  one  of  two  methods: 

•  By  prescribing  velocity,  (v11n),  or 
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Table  1 

MODEL  PARAMETERS  FOR  THE  PIVOTAL  EXPERIMENT 


Parameter 

Definition 

Value 

A 

eddy  viscosity 

150  m2  sec-1 

8 

(df/dy) 

2  X  10“11m"1sec"1 

fo 

Corollls  parameter 

8  X  10"5sec“1 

g' 

reduced  gravity  due  to 

stratification 

.02  m  sec"2 

Hi 

undisturbed  upper  layer  depth 

200  m 

h2 

undisturbed  lower  layer  depth 

reduced  gravity 

00 

2  layer 

2400  m 

Lx  X  Ly 

horizontal  dimensions  of  the 

model  domain 

420  km  x  280  km 

ax  X  y 

horizontal  grid  spacing  for  each 

dependent  variable 

10  km  x  5  km 

At 

time  step  reduced  gravity 

90Q  sec 

2  layer 

2700  sec 

v11n 

upper  leyer  Inflow  velocity 

30  cm/sec 

v21n 

lower  layer  Inflow  velocity 

.2  cm/sec 

a 

angle  of  Inflow 

21°  north  of  east 

*S 

Inflow  spin  up  time  constant 

30  days 

•  by  prescribing  transport,  (V^n),  and  allowing  the  model  to 
determine  the  velocities. 

Using  a  standard  Inflow  velocity  of  30  cm/sec  through  a  20  km  wide  200  m 
deep  port  determines  an  Inflow  transport  of  1.2  x  106  m/sec  or 
(106  m/sec  3  1  sv)  1.2  sverdrups.  This  value  Is  centered  within  the 
limits  of  the  observed  transport.  A  standard  Inflow  value  of  .2  cm/sec 
was  chosen  for  the  lower  layer  to  agree  with  observation.  Transport  out 
of  the  Medl terranean  Is  known  to  be  less  then  Incoming  transport  due  to 
the  effect  of  evaporation  being  greater  than  precipitation.  An  Inflow 
of  .2  cm/sec  through  a  245  km  wide  boundary  which  Is  a  maximum  of  2400  m 
deep  creates  a  mass  transport  of  1.1  sv.  It  also  creates  velocities 
along  the  southern  shelf  near  the  Strait  of  Gibraltar  of  10  cm/sec. 
This  also  agrees  with  observation  (personal  communication  -  J.C. 
Gascard).  Inflow  Is  spun  up  with  a  time  constant  of  thirty  days  to 
minimize  excitation  of  high  frequency  waves.  The  standard  ar  a  of 
Inflow  was  chosen  to  be  21  north  of  east  based  on  the  geometric 
orientation  of  the  Strait  of  Gibraltar.  This  angle  was  varied  based  on 
direct  observations  (Lacombe,  1961)  and  Inferences  from  satellite 
Imagery. 

The  wind  forcing  Is  derived  from  climatological  wind  stresses 
obtained  from  twenty  years  (1950-1970)  of  ship  observations  In  the 
Medl terranean  (May  1982).  Individual  stresses  were  estimated  from  ship 
observations  of  wind  speed  and  direction  using  a  quadratic  aerodynamic 
drag  law  with  a  drag  coefficient  dependent  on  the  wind  speed  and 
stability.  Monthly  averages  of  the  wind  stresses  were  then  calculated 
by  averaging  the  individual  wind  stress  estimates  from  each  month  on  a 
one  degree  latitude  by  one  degree  longitude  grid.  These  monthly 
averages  were  blllnearly  Interpolated  to  the  model  grid  (Figures  4,  5 
and  6).  Wind  forcing  In  the  model  Is  cyclic  with  a  period  of  one  year. 
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Monthly  mean  climatological  wind  stress  for  (a)  January,  (b)  February 
(d)  April. 
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Monthly  mean  climatological  Mind  stress  for  (a)  May,  (b)  June,  (c)  July  and 
(d)  August. 
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MODEL  RESULTS 


Model  results  will  be  presented  first  in  terms  of  two  layer 
results  then  explicit  reduced  gravity  results.  The  numerical  solutions 
presented  have  evolved  out  to  a  steady  state. 

The  pivotal  two  layer  experiment  used  the  parameters  of 
Table  1.  All  solutions  are  presented  in  terms  of  pynocllne  height 
anomaly  (PA)  the  deviation  of  the  Interface  between  the  upper  and  lower 
layers  from  Its  initial  rest  position  at  200  m  depth.  Solid  contours 
Indicate  a  deepening  of  the  layer  while  dashed  contours  Indicated  a 
shallowing  of  the  layer.  Closed  solid  contours  represent  regions  of 
antlcyclonlc  circulation  while  closed  dashed  contours  are  regions  of 
cyclonic  circulation  .  Figure  7a  shows  the  steady  state  solution  of  the 
pivotal  case.  The  dominant  feature  observed  Is  an  antlcyclonlc  gyre 
representing  the  Alboran  gyre  occupying  almost  the  entire  western  half 
of  the  basin.  This  solution  Is  quite  similar  to  the  semi -Implicit 
reduced  gravity  model  solutions  discussed  by  Preller  and  Hurlburt,  1982 
(Figure  7b).  The  addition  of  a  slowly  moving  lower  layer  seems  to  have 
very  little  effect  on  the  solutions.  Figure  8  shows  the  lower  layer 
velocity  field  for  the  pivotal  case.  Though  flow  enters  the  basin 
uniformly.  It  Is  quickly  steered  by  the  topography  to  the  southern  shelf 
and  then  exits  through  the  Strait  of  Gibraltar.  Although  the  velocity 
field  has  a  clearly  defined  pattern,  magnitudes  still  remain  very  small 
(<  10  cm/sec). 

A  more  Important  effect  seems  to  be  the  coastal  feature.  Cape 
Tres  Forces  and  the  Alboran  Island.  Figure  9a  shows  the  steady  state 
two  layer  solution  when  Inflow  Is  prescribed  due  east.  Cape  Tres  Forcas 
serves  to  deflect  the  jet,  not  allowing  It  to  abruptly  hit  the  southern 
boundary.  This  Is  brought  out  clearly  by  the  semi-implicit  reduced 
gravity  solution.  Figure  9b,  which  uses  rectangular  geometry.  The  jet. 
In  this  case,  hits  the  flat  wall  and  turns  sharply  upward  enhancing  the 
cyclonic  circulation  north  of  the  jet. 


LAYER  VEL0C 


A  series  of  experiments  was  also  made  to  test  the  effect  of  the 
Island.  Results  show  that  If  the  jet  entering  the  basin  at  an  angle  has 
too  large  a  velocity  (>  50  cm/sec),  the ‘Island  deflects  the  jet  to  its 
north  (Figure  10)  and  the  Alboran  gyre  disappears. 

A  number  of  experiments  were  also  performed  to  test  the  effect 
of  decreasing  Inflow.  Figure  11  shows  a  time  series  of  solutions  which 
was  Initialized  using  the  day  360  solution  of  the  pivotal  case. 
Velocity  was  then  allowed  to  decrease  to  1  cm/sec  over  a  period  of  180 
days.  Although  the  antlcyclonic  circulation  always  remains  In  the 
western  part  of  the  basin,  it  diminishes  both  in  strength  and  horizontal 
dimensions.  Also  of  note  is  the  small  cyclonic  cirulation  which  moves 
westward  and  intensifies.  This  movement  brings  about  the  existence  of 
two  gyres  west  of  Alboran  Island,  a  rather  unique  situation  sometimes 
observed  in  the  satellite  infrared  imagery. 

These  results  all  point  to  the  importance  of  knowing  the 
magnitude  of  the  Incoming  velocity  as  well  as  Its  direction.  The 
Alboran  gyre  appears  to  be  very  sensitive  to  these  changes.  Thus  In 
order  to  predict  the  characteristics  of  the  gyre,  a  time  series  of  the 
velocity  entering  the  Strait  of  Gibraltar  seems  necessary. 

Experiments  using  the  explicit  reduced  gravity  model  have  been 
recently  initiated  and  used  to  test  the  effects  of  climatological  wind 
stress.  This  test  Is  designed  to  determine  the  circulation  due  solely 
to  wind  stress  In  the  Alboran  Sea.  The  second  case  Is  forced  by 
30  cm/sec  Inflow  alone  for  360  days  and  wind  stress  force  Is 
then  added. 

Figure  12a,  b,  c  and  d  show  solutions  from  the  third  year  of 
Integration  at  mid  month  for  January,  April,  July  and  October.  The 
conclusion  may  be  drawn  from  these  solutions,  that  climatological  wind 
stress  alone  Is  responsible  for  a  somewhat  weak  but  consistent 
antlcyclonic  circulation  In  the  western  Alboran  Sea.  These  monthly 
solutions  also  show  that  wind  forcing  results  in  small  time  variations 
In  the  circulation  of  the  Alboran  Sea.  Figure  13  represents  the  basin 
averaged  kinetic  and  potential  energy  for  the  three  years  of 
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Figure  13.  Basin  averaged  kinetic  and  potential  energy  for  the  three 
years  of  the  wind  driven  case. 


Integration.  These  curves  show  that  although  energy  varies  during  the 
year,  a  yearly  pattern  develops  and  repeats.  Indicating  that  a  steady 
state  has  been  reached. 

Figure  14  presents  the  solution  of  the  second  experiment  after 
350  days  of  Integration  using  30  cm/sec  Inflow.  At  day  360,  the 
climatological  wind  forcing  was  added  and  the  model  was  Integrated  for 
three  more  years.  Results  from  January,  April,  July  and  October  of  the 
fourth  year  are  presented  In  Figure  15a,  b,  c  and  d.  The  addition  of 
the  wind  enhances  the  horizontal  dimensions  and  Intensity  of  the  gyre. 
These  solutions  also  show  small  scale  changes  In  the  circulation  with 
time  due  to  the  wind  stress.  Figure  16  shows  the  basin  averaged  energy 
curves  for  the  four  years  of  this  experiment.  Solutions  appear  to 
approach  a  steady  state  by  year  four.  Thus  from  a  climatological  point 
of  view  wind  stress  assists  the  development  of  an  anti cyclonic  gyre 
throughout  the  year. 
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INTERFACE  DEV  I  ATI  0N  alb  stress 

DH=  10  M  DAY=  350 


Solution  of  the  port  driven  (30  cm/sec)  only,  explicit  reduced  gravity  case 


INTERFACE  DEV1RT10N  alb  stress  INTERFACE  DEV  I  ATI  0N  alb  stress 


Figure  15.  Solutions  of  the  wind  and  port  driven  explicit  reduced  gravity  case  for 
(a)  January,  (b)  April,  (c)  July  and  (d)  October  of  the  fourth  year. 
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Figure  16.  Basin  averaged  kinetic  and  potential  energy  for  the  four 
years  of  the  port  and  wind  driven  case. 


SUMMARY 


The  non-linear  reduced  gravity  and  two  layer  models  of  Hurlburt 
and  Thompson  (1980)  have  been  adapted  (Preller  and  Heburn,  1983)  to 
study  circulation  In  the  Alboran  Sea.  Model  results  using  real  geometry 
and  topography  In  a  two  layer  model  show  little  effect  upon  the  upper 
layer  by  the  slowly  moving  lower  layer.  More  Important  are  the  effects 
of  a  coastal  feature  Cape  Tres  Forcas,  which  deflects  the  flow  from  the 
southern  boundary.  Also,  the  Alboran  Island  deflects  the  Incoming  jet 
to  Its  north  If  velocities  greater  than  50  cm/sec  are  used  at  Inflow. 
Additional  experiments  show  the  high  sensitivity  of  the  Alboran  gyre  to 
changes  In  the  magnitudes  of  the  Inflow,  Indicating  that  a  thorough 
knowledge  of  the  Inflow  magnitude,  direction  and  time  variation  Is 
necessary  to  predict  the  characteristics  of  the  gyre.  Reduced  gravity 
solutions  forced  by  monthly  mean  climatological  wind  stress  shows  that 
climatological  wind  stress  enhances  the  gyre  throughout  the  year  and 
causes  some  small  time  varying  changes  In  the  circulation  of  the  entire 


REFERENCES 


Bethoux,  J.  P.,  1979.  Budgets  of  the  Medl terranean  Sea.  Their 

dependence  on  the  local  climate  and  on  characteristic  of  the 
Atlantic  waters.  Oceanol.  Acta.  2,  2:  137-163. 

Burkov,  V.  A.,  Drivosheya,  V.  G.  Ovchinnikov,  I.  M.  and  Savin,  M.  T., 
1979.  Eddies  In  the  current  system  of  the  western 

Medl terranean  Basin.  Oceanology,  19:  9-13. 

Cheney,  R.  E.,  1977.  Recent  observations  of  the  Alboran  Sea  front. 
NAVOCEANO  Technical  Note  370073-77,  24  pp. 

Cheney,  R.  E.,  1978.  Recent  observations  of  the  Alboran  Sea  frontal 
system.  J.  Geophys.  Res.,  83:  4593-4597. 

Gallagher,  J.  J.,  Fecher,  M.  and  Gorman,  J.,  1981.  Project  HUELVA 
oceanographic/acoustic  Investigation  of  the  western  Alboran 
Sea.  NUSC  Technical  Report  6023,  106  pp. 

Grousson,  R.  and  Faroux,  J.  1963.  Measure  de  courants  de  surface  en 
Mer  d' Alboran.  Cah.  Oceanogr.,  15:  716-721. 

Hurlburt,  H.  E.  and  Thompson,  J.  D.,  1980.  A  numerical  study  of  Loop 
Current  Intrusions  and  eddy  shedding.  J.  Phys.  Oceanogr.,  10: 
1611-1651. 

Katz,  E.  J.,  1972.  The  Levantine  Intermediate  water  between  the 
Straight  of  Sicily  and  the  Strait  of  Glbralter.  Deep  Sea  Res., 
19:  507-520. 

Lacombe,  H.,  1961.  Contribution  A  L'Etude  du  Regime  du  Detroit  de 

Glbralter.  Cah.  Oceanogr.,  XIII:  74-107. 

Lacombe,  H.,  1971.  Le  Detroit  de  Gibraltar.  Note  et  Memolres  de 

Service  Geologlque  du  Maroc.,  No  222  bis,  111-146. 

Lacombe,  H.,  1982.  Regime  of  the  Strait  of  Gibraltar  and  of  Its  east 
and  west  approaches.  In:  J.  C.  J.  Nlhoul  (Editor), 

Hydrodymanlcs  of  semi-enclosed  seas.  Elsevier,  Amsterdam,  pp. 

Lanolx,  F.,  1974.  Project  Alboran  Etude  Hydrologlque  Dynamlque  de  la 
Mer  d'Alboran.  Tech,  report  66,  N.  Atl .  Treaty  Org.  Bussels, 
pp  39. 

May,  P.W.,  1982.  Climatological  flux  estimates  In  the  Mediterranean 
Sea:  Part  I.  Winds  and  wind  stress.  Naval  Ocean  Research  and 
Development  Activity,  NSTL  Station,  MS.  NORDA  Technical  Report 
54,  p.  56. 


Mommsen,  D.  B.,  Jr.,  1978.  The  effect  of  wind  on  sea  surface 
temperature  gradients  in  the  Strait  of  Gibraltar  and  Alboran 
Sea.  Report  from  Fleet  Weather  Central,  Rota,  Spain,  18  pp. 

Nof,  0.,  1978.  On  geostrophlc  adjustment  in  Sea  straits  and 

estuaries:  theory  and  laboratory  experiments.  Part  II.  Two- 
layer  system  J.  Phys.  Oceanogr.,  8:  861-872. 

Ovchinnikov,  I.  M.,  1966.  Circulation  in  the  surface  and  intermediate 
layers  of  the  Mediterranean.  Oceanology,  6:  48-59. 

Ovchinnikov,  I.  M.,  Krlvosheya,  V.  G.  and  Maskalenko,  L.  V.,  1976. 

Anamalous  features  of  the  water  circulation  of  the  Alboran  Sea 
during  the  summer  of  1962.  Oceanology,  15:  31-35. 

Peluchon,  G.  and  Oonguy,  J.  R.,  1962.  Travauz  Oceanographiques  d 
"l'Orlgny"  dans  le  Detroit  de  Gibraltar.  Compalgne 
Internationale  -  15  mpl,  15  Juln  1961.  Zeml  partle. 

Hydrologle  en  Mer  d'Alboran.  Cah.  Oceanogr.,  14:  573-578. 

Petit,  M.,  Klaus,  V.  Gelcl,  R.,  Fusey,  F.,  Thery,  J.  J.  and  Bouly,  P., 
1978.  Etude  d'un  tourblllon  oceanlque  d'echelle  moyenne  en  mer 
d'Alboran  par  emplol  conjoint  techniques  spatlales  et 

oceangraphlques.  C.  R.  Acad.  Scl.  287:  215-218. 

Porter,  D.  L.,  1976.  The  antlcyclonlc  gyre  of  the  Alboran  Sea, 

Independent  research  report  from  M.I.T.-WH0I  joint  program. 
Woods  Hole,  Ma,  pp  29. 

Stevenson,  R.  W.,  1977.  Huelva  Front  and  Malaga,  Spain  eddy  chains  as 
defined  by  satellite  and  oceanographic  data.  Deut.  Hydrogr.  A. 
30,  2:  51-53. 

Stommel,  H. ,  Bryder,  H.  and  Magelsdorf,  p.,  1973.  Does  some  of  the 
Mediterranean  outflow  come  from  great  depth?  Pure  and  Appl 
Geophysics.  105:  879-889. 

Wannamaker,  B.,  1979.  the  Alboran  Sea  Gyre:  Ship  satellite  and 

historical  data  SCAIANT  ASW  Research  Centre  Report  Sr-30,  La 
Spezla,  Italy,  27  pp. 

Preller,  R.  and  H.  E.  Hurlburt,  1982.  A  Reduced  gravity  numerical  model 
of  circulation  In  the  Alboran  Sea.  In:  J.C.J.  Nihoul  (Editor) 
Hydrodynamics  of  semi -enclosed  seas.  Elsevier,  Amsterdam,  75- 
89. 


Preller,  R.  H.  and  G.  W.  Heburn,  1983.  A  simple  numerical  model  of 
circulation  In  the  western  Mediterranean  Sea  submitted  to 
Oceanollglca  Acta. 


Vanwyckhouse,  R.  J.,  1973.  Synthetic  Bathymetric  Profiling  System 
(SYNBAPS).  Naval  Oceanographic  Office,  Washington,  DC 
Technical  Report  233,  p.  138. 

Vanwyckhouse,  R.  J.,  1979.  SYNBAPS,  Volume  I  -  Data  sources  and  data 
preparation.  Naval  Ocean  Research  and  Development  Activity, 
NSTL  Station,  Ms,  NORDA  Technical  Note  35. 


VERTICAL  SHEAR  FROM  TOPS 


INTRODUCTION 

The  growth  and  propagation  of  Internal  waves  In  the  world 

oceans  Is  often  attributed  to  the  regions  of  strong  gradients  In 

velocity  (vertical  shear)  and  density.  To  predict  the  spatial  and 
temporal  distribution  of  density  and  velocity  gradients  would  require  an 
extensive  effort  In  ocean  surveys  modeling. 

Numerous  field  experiments  have  used  high  resolution  velocity 
and  velocity/temperature  microstructure  profilers  to  measure  vertical 
variation  of  shear,  temperature,  Richardson  number  and  turbulence 
dlsspatlon.  For  a  review  of  these  experiments  see  Harding  et  al., 
1983.  Though  these  experiments  were  performed  at  various  regions  in  the 
northern  hemisphere,  most  were  performed  over  relatively  short  time 
scales  (a  few  days).  Thus  It  would  be  a  practical  impossibility  to 

measure  or  predict  from  experimental  data  sources  alone,  the  density  and 

velocity  gradients  over  large  time  and  horizontal  space  scales. 
Numerical  forecasting  of  these  upper  ocean  properties  with  verification 
by  properly  designed  field  programs  serves  as  a  useful  alternate 
approach. 

The  best  candidate  for  such  a  modeling  effort  Is  the 
Thermodynamic  Ocean  Prediction  System  (TOPS)  currently  operational  at 
Fleet  Numerical  Ocean  Center  (FNOC)  (Clancy  and  Martin,  1981).  This 
model  provides  the  upper  ocean  currents  required  to  synoptlcally 
forecast  upper  ocean  vertical  shear  over  large  regions  of  the  world's 
ocean.  TOPS  additionally  provides  velocity  components  through  the  mixed 
layer  making  It  potentially  useful  for  shear  forcastlng. 


MODELING  APPROACH 


The  TOPS  model  Is  basically  a  Level  II  turbulence  closure  model 
of  Mellor  and  Yamada  (1974).  The  model  results  In  an  isopycnal  mixed 
lAyer  with  vertical  uniformity  due  to  mixing  caused  by  wind  Induced 
shear  Instability  and  surface  cooling.  During  strong  wind  forcing 
and/or  strong  surface  cooling,  the  mixed  layer  deepens  due  to 
entrainment  from  below  caused  by  turbulent  mixing.  During  weak  wind 
forcing  and/or  strong  surface  heating,  the  turbulent  energy  decreases  to 
the  point  where  active  entrainment  at  the  base  of  the  mixed  layer  Is  no 
longer  maintained  and  the  layer  shallows.  This  model  has  been 
extensively  tested  against  observations  by  Warn-Varnas,  Dawson  (1981) 
and  Martin  (1982,  1983)  at  specific  experimental  and  weather  ship 
locations.  It  has  also  been  tested  In  the  operational  sense  by  Clancy 
and  Martin  (1979),  Clancy  et  al.,  (1981),  and  Clancy  and  Poliak  (1982). 

Forcing  for  the  TOPS  model  Is  derived  from  the  Navy's  global 
weather  prediction  model  NOGAPS  In  the  form  of  surface  stresses  and  heat 
fluxes.  It  Is  Initialized  and  updated  by  an  analysis  scheme  which  uses 
the  TOPS  forecast  products  for  first  guess  fields  and  blends  them  with 
dally  XBT  observations  and  satellite  derived  SST  fields.  The  TOPS  model 
Is  essentially  a  3-D  mixed  layer  model  consisting  of  a  grid  of  profiles, 
but  with  advectlon  of  temperature  and  salinity  between  grid  points.  At 
FNOC,  TOPS  runs  on  the  standard  63  x  63  northern  Hemispheric  polar 
stereographic  grid  with  381  km  grid  resolution  at  60  N,  Figure  17.  The 
vertical  grid  has  Irregularly  spaced  levels  with  az  varying  from  5  km  at 
the  surface  to  100  km  below  the  200  m  level  (Figure  18).  Bottom 
boundary  conditions  specify  constant  values  of  temperature,  salinity  and 
velocity;  surface  boundary  conditions  specify  fluxes  of  heat, 
precipitation  or  evaporation  and  momentum.  Lateral  boundary  conditions 
specify  no  Inflow  or  outflow  and  zero  fluxes  of  heat  and  salinity. 
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Figure  18.  Vertical  grid  utilized  by  TOPS.  The  quantities  T,  S,  u,  v, 
u.  and  v  are  defined  at  the  depths  indicated  in  the  figure. 

a  a 

All  turbulence  quantities  and  w.  are  defined  at  depths  midway 

a 

between  those  shown  in  the  figure. 


EXPERIMENTAL  APPROACH 


The  purpose  of  this  study  was  to  continue  the  research 
Initiated  last  year  (Harding  et  al.,  1983)  to  investigate  the  potential 
of  TOPS  for  shear  prediction  and  to  look  at  temporal  and  spatial  vari¬ 
ability  of  the  TOPS  output.  In  particular,  this  study  Investigated  long 
term  shear  variations  as  well  as  seasonal  comparisons. 

At  the  time  of  these  experiments,  velocity  fields  from  TOPS  had 
not  been  archived  for  a  long  enough  period  of  time  to  examine  long  term 
variations.  Thus  the  version  of  TOPS  which  resides  on  the  Texas  Instru¬ 
ment  Advanced  Scientific  Computer  (TIASC)  at  the  Naval  Research  Labora¬ 
tory  (NRL)  In  Washington  was  run  continuously  In  a  hlndcast  mode  for  six 
months.  This  version  of  the  model  is  described  In  detail  by  Warn-Varnas 
et  al.,  (1983)  who  used  It  to  predict  large  scale  thermal  anomalies  in 
the  North  Pacific.  The  FNOC  analysis  and  model  outputs  required  to 
force  TOPS  are  also  available  at  NRL  over  the  entire  northern  hemisphere 
for  extended  time  periods. 

The  chosen  area  of  study  Is  12  x  21  point  subregion  of  the  63  x 
63  FNOC  grid  over  the  North  Pacific,  Figure  19.  The  TOPS  model  was 
Integrated  over  the  six  month  period  March  2  through  August  29,  1981. 
The  model  Is  Initialized  with  March  climatological  salinity  and  March  2 
thermal  structure.  The  model  was  then  run  for  approximately  six  months 
using  FNOC  6  hourly  wind  and  heat  flux  fields.  A  detailed  explanation 
of  initial  conditions,  boundary  conditions  and  forcing  functions  are 
found  In  Warn-Varnas  et  al.,  (1983).  An  experiment  has  been  conducted 
to  determine  over  what  period  of  time  the  solutions  are  affected  by 
spurious  motions  Induced  by  Initialization.  This  time  period  appears  to 
be  15-20  days.  Thus  only  solutions  after  a  15  day  Initialization  period 
are  considered.  The  FNOC  Extended  Ocean  Thermal  Structure  (EOTS)  pro¬ 
duct  provided  the  thermal  structure  Initialization  fields  required.  Ta¬ 
ble  2  lists  the  forcing  fields,  their  source  and  frequency. 

From  hourly  TOPS  output,  scalar  shears 


S  3  ((du/dz)2  +  (dv/dz)2f/2 


were  calculated  at  levels  of  TOPS  located  velocities.  Z  represents  the 
vertical  distance  between  the  velocities  given  by  the  model.  Also 
calculated,  was  the  gradient  Richardson  number, 

R1  =  -(gdP/dz)/S2 

( p  3  water  density  at  a  given  level,  g  s  gravitation  acceleration)  at 
levels  Intermediate  to  the  velocity  levels.  This  number  is  a  measure  of 
the  stability  of  the  water  column  at  that  level.  Yamada  (1975),  using  a 
turbulence  closure  formulation  similar  to  that  used  in  TOPS,  determined 
the  critical  flux  of  Richardson  number  exists  between  .18  -  .27.  Values 
above  .25  indicate  stability  In  the  water  column.  Thus  for  this  study 
values  between  0.  and  .25  will  be  used  to  indicate  shear  instability. 


FNOC 

TABLE  2 

FIELDS  USED  FOR  TOPS  FORCING 

FIELD 

SOURCE 

FREQUENCY 

Boundary  Layer  Wind 
Speed 

Planetary  Boundary 

Layer  Model 

6  hourly 

Boundary  Layer 

Wind  Direction 

Planetary  Boundary 

Layer  Model 

6  hourly 

Solar  radiation  Flux 

Primitive  Equation 

Model 

6  hourly 

Total  Heat  Flux 

Primitive  Equation 

Model 

6  hourly 

Evaporative  Heat  Flux 

Primitive  Equation 

Model 

6  hourly 

Precipitation 

Primitive  Equation 

Model 

12  hourly 
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TOPS  model  results  will  be  presented  in  terms  of  contour  plots 
and  frequency  distributions  of  shear  and  Richardson  number  for  various 
levels  and  vertical  cross  sections  within  the  North  Pacific  domain.  The 
purpose  Is  to  investigate  spacial  and  temporal  variability  of  the 
potential  FNOC  shear  product.  The  vertical  cross  sections  of  Richardson 
number  show  the  dividing  line  between  stable  and  unstable  water,  the 
pycnocline.  Regions  of  shear  Instability  are  shaded  in  these  figures. 
Note  that  In  all  plots  depicting  shear,  values  have  been  scaled  upward 
by  a  factor  of  104. 


RESULTS 


In  previous  studies  (Harding  et  al.,  1983),  we  examined  the 
Importance  of  short  time  scales,  six  hourly  variations,  diurnal 
variations  and  dally  variations.  These  results  showed  a  high  degree  of 
variability  even  on  the  shortest  time  scales.  In  this  study,  we  will  be 
concerned  with  seasonal  time  scales  and  long  term  trends  in  shear. 

The  first  series  of  contoured  plots  will  show  monthly  values  of 
shear  for  the  six  months  studied.  Values  are  contoured  for  various 
levels  and  cross  sections  starting  on  March  17  and  then  depicted 
approxlmatley  every  30  days. 

Figure  20  shows  lines  of  latitude  and  longitude  along  with  land 
points  over  the  North  Pacific  region  studied.  All  level  plots  will  be 
presented  with  this  orientation.  With  reference  to  this  horizontal 
orientation,  cross  sectional  plots  are  made  along  lines  drawn  from  the 
center  point  at  the  left  boundary  (6  of  12)  and  the  center  point  of  the 
top  boundary  (12  of  24). 

Figure  21  shows  a  monthly  time  series  of  shear  at  level  2  (7.5 
m).  Most  noticeable  In  this  series  of  figures  Is  the  fact  that  the 
shear  Intensifies  from  March  to  August.  The  upper  left  corner  of  these 
plots  develops  Into  a  region  of  consistently  high  shear  (200  x  10“4). 
Geographically,  this  region  Is  to  the  south-southwest  of  the  Bering 
Sea.  It  Is  a  region  of  frequent  Intense  storm  and  frontal  passage. 
Harding  et  al.,  1983,  have  shown  that  shear  Instability  appears  directly 
correlated  to  such  frontal  passage. 

Figures  22  and  23  are  histograms  of  the  shear  at  level  2  (7.5 
m)  for  the  corresponding  times.  These  figures  clearly  show  the  Increase 
In  magnitude  of  shear  with  time.  Initially  the  average  shear  value 
centers  around  8  x  10”3.  By  summer,  the  average  value  is  1.4  x  10"2 
with  much  higher  values  appearing  then  in  spring. 

The  same  times  series  for  the  long  (left  to  right)  cross 
section  Is  shown  In  Figure  24.  This  time  series  shows  the  same  trend  as 
In  the  level  series,  larger  magnitude  shear  developing  In  the  summer. 


Figure  21.  Monthly  variation  of  the  7.5  m  shear  (s"1  x  10 
months  March  through  August. 
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Figure  23.  Frequency  distribution  of  the  7.5  m  shear  (s"1  x  IQ4)  for 
the  months  June  through  August. 
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Figure  24.  “Long"  monthly  cross  sections  of  shear  (s”1  x  104)  for  March 
through  August. 
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However,  It  shows  one  more  very  important  trend.  If  spring  is  defined 
as  March-Aprll  with  May -June  as  the  transition  period  into  summer,  July- 
August,  it  becomes  clear  that  shear  extends  to  much  greater  depth  in 
spring  then  in  summer.  The  deepest  shear  lies  in  the  region  south  of 
the  Bering  Sea,  the  same  trend  which  was  observed  in  the  horizontal 
sections.  The  shallowing  tendency  is  also  apparent  in  the  shorter  cross 
sections.  Figure  25.  Again  shear  Intensifies  at  the  surface  in  the 
summer  but  extends  to  much  greater  depths  in  spring.  Figures  26  and  27 
are  histograms  of  shear  at  level  six.  These  figures  also  show  the 
dramatic  decrease  in  shear  from  spring  to  summer  at  32.5  m  depth.  (Note 
the  difference  in  scale  on  the  y-axls  when  comparing).  This  existence 
of  shear  at  depth  In  late  winter  -  early  spring  Is  not  an  unexpected 
phenomenon.  Figure  19  shows  the  location  of  Ocean  Station  Papa  and 
Ocean  Station  November  while  Figures  28  and  29  show  values  of  mixed 
lawyer  depth  during  1961  obtained  from  these  stations.  The  maximum  mixed 
layer  depth  exists  In  March  and  slowly  decreases  until  there  Is  a 
dramatic  shallowing  In  June,  July  and  August.  This  well  known  spring 
transition  of  the  mixed  layer  occurs  abruptly  In  response  to  an  extended 
period  of  weak  winds  and  strong  surface  cooling.  Figure  30  further 
demonstrates  the  decrease  In  shear  at  depth  using  horizontal  sections  of 
the  model  results  at  level  6  (32.5  m).  Active  shear  exists  at  this 
depth  In  March  and  April,  diminishing  In  May  and  all  but  disappearing 
Into  the  summer  months.  This  shallowing  of  the  mixed  layer  also 
explains  the  Increase  In  shear  magnitude  during  the  summer.  The  model 
formulation  conserves  momentum  within  the  mixed  layer,  thus  as  the  layer 
shallows,  velocity  as  well  as  shear  will  Increase  within  the  narrow 
mixed  layer. 

Level  6  was  chosen  as  an  average  depth  at  which  shear  might 
exist  for  the  entire  six  month  period.  The  maximum  depth  to  which  shear 
extends  Is  shown  clearly  In  the  vertical  cross  sections.  Figure  31 
shows  values  of  shear  starting  at  noon  of  March  17  and  then  for  forty- 
eight  hour  Intervals  until  noon  of  March  23.  Using  a  contour  Interval 
of  5,  the  maximum  depth  of  these  contours  extend  below  125  m.  Figure  32 
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Figure  25.  "Short"  monthly  cross  sections  of  shear  (s_1  x  104 
March  through  August. 
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Figure  26.  Frequency  distribution  of  32.5  m  shear  (s-1  x  104)  for  March 
through  May. 
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represents  Richardson  numbers  for  this  same  time  period.  These  show  the 
region  of  Instability  extending  down  to  approximately  100  m  depth. 
These  two  figures  show  how  closely  regions  of  Instability  correspond  to 
regions  of  shear.  Figure  33  shows  the  same  time  series  of  shear  for 
April  (16-22).  Values  of  shear  extend  to  approxlmatley  60  m  depth, 
having  shallowed  somewhat  from  their  March  position.  Richardson  number 
cross  sections  (Figure  34)  for  this  period  show  the  development  of  an 
Interesting  feature,  two  vertically  separate  regions  of  Instability. 
The  shallower  region  extends  from  the  surface  to  25  m  depth  while  the 
deeper  region  exists  at  an  average  depth  of  80  m.  This  secondary  deep 
Instability  Is  a  remnant  of  the  deeply  mixed  layer  which  existed  In 
March. 

In  our  previous  study  of  shear,  a  strong  correlation  existed 
between  frontal  passage  and  dynamic  stability  of  the  mixed  layer. 
Hal  pern  (1974)  saw  the  same  correlation  from  moorings  in  the 
northeastern  Pacific.  Figures  35  and  36  are  examples  of  fronts  overlain 
on  Richardson  number  plots  calculated  at  10  m  depth  In  March  and 
April.  Most  of  the  area  In  these  figures  represents  a  Richardson  number 
between  0-.25.  Regions  of  frontal  passage  are  related  to  Richardson 
number  values  between  0-.25  while  values  greater  than  .25  are  found  In 
regions  of  high  pressure  systems. 

Figures  37  and  38  show  histograms  of  Richardson  number  at  10  m 
for  the  thirty  day  Intervals  presented  In  Figure  20.  Once  again,  note 
the  differences  of  scales  on  the  y-axis.  The  percentage  of  Richardson 
number  lying  between  0-.25  drops  dramatically  In  June,  July  and  August. 

A  seasonal  comparison  of  shear  Is  depicted  In  the  following 
figures.  A  sixty  day  TOPS  model  run.  Initialized  at  October  29,  1980, 
was  conducted  for  November  and  December  of  1980.  Solutions  at  noon  of 
November  13,  15,  17  and  19  are  presented  in  comparison  with  solutions  at 
noon  of  May  16,  18,  20  and  22.  Horizontal  cross  sections  at  7.5  m. 
Figures  39  and  40,  for  these  two  time  periods  show  that  shear  near  the 
surface  Is  higher  during  the  spring  then  during  the  fall.  In  the  upper 
left  hand  corner  of  these  plots,  that  region  south-southwest  for  the 
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Figure  34.  "Short"  cross  sections  of  Richardson  number  for  the  same 
time  period  as  that  of  Figure  33. 


Figure  35.  Level  2  (7.5  m)  Richardson  number  plots  with  atmospheric 
fronts  overlaid  for  the  same  time  period  as  Figure  31. 
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2  (7.5  m)  Richardson  number  plots  with  atmospheric 
overlaid  for  the  same  time  period  as  Figure  33. 
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Figure  37.  Monthly  Level  2  (7.5  m)  Richardson  number  frequency 

distribution  for  March  through  May. 
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Figure  38.  Monthly  Level  2  (7.5  m)  Richardson  number  frequency 

distribution  for  June  through  August. 
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beginning  at  noon  on  May  16 


Bering  Sea  values  of  shear  are  twice  as  large  in  May  as  in  November. 
The  vertical  cross  sections,  Figures  41,  42,  43  and  44  also  indicate 
that  shear  is  higher  in  spring  than  fall,  and  deeper  in  fall  than 
spring.  This  comparison  agrees  with  the  mixed  layer  depths  at  stations 
Papa  and  November.  Both  of  these  stations  indicated  that  the  mixed 
layer  depth  Is  greater  In  November  than  in  May.  Solutions  at  an  average 
depth  of  32.5  m  (level  6)  for  both  seasons.  Figures  45  and  46,  compare 
quite  closely. 
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"Short"  cross  sections  of  shear  (s-1  x  104)  for  the  48  hour 
Intervals  beginning  at  noon  on  May  16  through  noon  on  May 
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“Long"  cross  section  of  shear  (s"1  x  1(F)  for  the  48  hour 
Intervals  beginning  at  noon  on  November  13  through  noon  on 
November  19. 
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SUMMARY  AND  CONCLUSIONS 


N. 


This  study  has  attempted  to  simulate  the  statistics  of  temporal 
and  geographical  distribution  of  shear  over  the  upper  400  m  of  the 
world's  oceans.  In  particular  we  have  examined  a  large  area  of  the 
North  Pacific  over  the  six  month  period  March  through  August  of  1981. 
Forcing  was  provided  by  the  prediction  fields  obtained  from  FNOC's 
hemispheric  weather  model. 

Model  resul ts  have  shown : , 


•  "Surface  shear  Increases  from  late  winter  -  early  spring  to 
^ummer\^ 

•  The  region  south-southwest  of  the  Bering  Sea  develops  Into  a 

region  of  very  high  surface  shear  In  spring  and  summer. 
a  — . -  •  . .  .  j 

•  The  depth  to  which  shear  extends  Is  closely  correlated  to 

the  mixed  layer  depth,  thus  In  winter  shear  values  are  found 

at  greater  depths  then  In  summer. 

s.  ■ ■  _  >  > 

•  Seasonal  comparisons  show  that  the  magnitude  of  shear  Is 

larger  In  May  then  In  November,  but  extends  to  greater 
depths  In  November  then  In  May. 

•  ^hear  Instability  Is  strorfgly  linked  to  the  passage  of 

atmospheric  fronts.  Regions  of  stability  occur  under  high 
pressure  systems  and  low  winds. 

<4  . . . . . *  . .  t  - 

•Very  little  convective  instability  exists  during  the  six 
month  period  examined  (Rpn  <  0). 
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